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The problem of propagation of vertically polarized radiowaves in an 
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that the Fourier kernels of the transform equations in the split step algorithm represent 
planes waves and (ii) compensating for the effects of rough ocean surfaces by using a rough 
surface reduction factor directly in the spectral domain. To accomplish this a redefinition of 
the Fourier transform pair is done to ensure mathematical consistency. The formulation also 
incorporates the first and second derivatives of the reffactivity index to accommodate steep 
gradients in the reffactivity profile. Hanning windows are used in both the spatial and 
wavenumber domains to contain computational requirements. The effects on propagation 
by varying parameters such as wave heights, computational domain ceilings, frequency and 
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I. INTRODUCTION 



A. BACKGROUND 

Propagation of radiowaves over rough ocean surfaces and in an inhomogeneous 
atmosphere is a topic of particular interest to the Navy. The effects of this type of 
environment on the strength of a signal at a given distance downrange from a transmitter are 
significant, and play a critical role in determining if a communication link can be 
successfully established, or a target can be detected by radar. This type of radiowave 
propagation is governed by a modified Helmholtz equation which is an elliptic partial 
differential equation. Solving this type of equation is computationally very intensive and 
becomes impractical when addressing typical propagation problems. The Helmholtz 
equation can be approximated by a parabolic equation (PE) if only one way propagation is 
considered, which is a reasonable assumption for most typical propagation problems. The 
advantage of using the parabolic equation approach is that a solution can be obtained much 
more efficiently through the use of a range stepping technique. Kuttler and Dockery [3] 
discuss the basic idea and various approximations involved in the development of the PE. 
This technique makes it possible to easily estimate propagation losses several hundreds of 
kilometers downrange with antenna heights up to several hundreds of meters for 
frequencies through the Super High Frequency (SHF) band. By only considering waves 
propagating in the forward direction, the strength of a signal at a given location downrange 
will be determined by direct and reflected (from the ocean surface) waves. Interference 
between these two waves will result in reflection multipath fading which can result in high 
gain or severe loss of signal. The effects of reflection multipath fading will depend greatly 
upon the surface height deviations of the ocean surface which will be determined by wind 
speed. 
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B. OBJECTIVE 



In this thesis the split-step PE algorithm is used to predict propagation of vertically 
polarized radiowaves over rough ocean surfaces in the presence of an inhomogeneous 
atmosphere. This thesis builds upon work previously done by Janaswamy [1] in this area 
which addressed horizontal polarization. In this approach the Fourier transform pair used in 
the split-step PE algorithm is modified to accommodate the effects of sea surface roughness 
directly into the formulation. The idea behind this approach will be to cast the transform 
equations in terms of incident and reflected waves, and then use the rough surface reduction 
factor available for plane waves, according to Miller [2], directly in the spectral domain. 
Chapter II presents the derivation and formulation of the modified Fourier transform pair 
for the split-step PE algorithm. Chapter III details the generation of the numerical procedure 
for solving the parabolic equation. In Chapter IV the performance of the numerical solution 
is examined by varying several important parameters and observing the effects on 
radiowave propagation. The effects of changing parameters such as step size, atmospheric 
refractivity profile, number of points for Fast Fourier Transform (FFTs), frequency, and 
wind speed will be investigated. Recommendations and conclusions are presented in 
Chapter V. 
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II. FORMULATION 



In this chapter we present the theory governing the parabolic equation and the 
derivation of the Fourier transform pair to be used with the split-step algorithm. This 
derivation parallels that given by Janaswamy [1] which considered the same propagation 
problem discussed here, but for horizontally polarized radiowaves. This chapter reiterates 
many of the formulas from that report but modifies them as appropriate to derive correct 
formulas for the case of vertically polarized waves. 

Figure 1 illustrates the basic geometry of the problem we are investigating. We 
consider the source at an initial range (x = 0) and height (z = z,) to be an omnidirectional 
point source. Given the position of the receiver, frequency of operation, refractive index 
profile of the atmosphere, wind speed, and the ground constants (s r and a), we wish to 
determine the signal strength along the path between the transmitter and receiver. 

^max 



r 




Figure 1. Source Producing Fields Over Rough Ocean Surface 
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To solve this propagation problem the standard parabolic equation as defined below 
by Kuttler and Dockery [3] is used: 



d 2 u\x,z) , duix.z) , r / \ 1 / \ 

— + 2 ik 0 — ^ + 2k~ \m(x, z) - l]«(x, z) = 0 , 



( 1 ) 



where u{x,z) = refers to a spherical coordinate reference system. 

fl.0%9) is the ^-component of the magnetic field. Equation (1) assumes an e '“'time 
dependence. The coordinate system used in this thesis is shown in Figure 2 below. 



z 




Figure 2. Earth Centered Spherical Geometry 
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In this coordinate system the transmitter is located at 9 = 0 and r = a e + z, where a t is the 
radius of the earth and z t is the height of the transmitter above the surface of the earth. The 

free space wave number is given by k 0 = oo^/e 0 p 0 ,x is the range axis, z is the height axis 
and n^x, z) is the modified refractive index under earth flattened conditions as given by 

Kuttler and Dockery [3] and is equal to: m(x,z) = , with n being the actual refractive 

index. 

By examining (1) we can see that the highest derivative in terms of the range, x, is 
of the first order for the parabolic case, unlike a Helmholtz equation which would be of 
second order. This desirable characteristic of the parabolic equation allows us to use (1) to 
solve for a field at a given range based upon a known field at a previous range. If the range 
step size, or Ax, is kept reasonably small and the refractive index varies slowly, the field at a 
new range can be computed from the field at a previous range by using the split-step 
algorithm by Tappert [4] below: 





-ip 2 Ax 1 


u(x,z) = e ik ° {m -' ) * x F~ l < 


e 2k ° f[m(x 0 ,z)] 



Later in this thesis we will modify this equation to accommodate refractive indices which 
have steep gradients. F represents the forward Fourier transform operator corresponding to 
an appropriate spectral decomposition of the field in the vertical direction, and F ' 1 is its 
inverse. 

To develop the form for the Fourier operators to be used in the split-step algorithm 
for rough ocean surfaces it is logical to first solve for the smooth surface case, and then 
modify this result so it can be applied to a rough ocean surface. In our formulation the ocean 
surface shows up as a boundary condition for which the parabolic equation must be solved. 
For the smooth ocean surface case the parabolic equation as given at (1) must be solved 
subject to the impedance boundary: 
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3mx,0) , x 

— - + a 0 u(x,0) = 0 , 



where a 0 = 



ik oy [e~-\ 



. The complex dielectric constant, s re is defined as: 



( 3 ) 



£ = £ + lO . 

rc r r ? 



where, cr r is the relative conductivity of the earth given by: 

cr 



<*r = 



CO£r 



to = angular frequency (rads / s) 



(4) 



(5) 



The boundary condition given at (3) is valid only for |s rc | » 1 , which for ocean surfaces is 
easily met. To solve the standard parabolic equation given at (1) subject to the boundary 
condition at (4), the following mixed Fourier transform pair by Kuttler and Dockery is used 

[3]: 

oo 

u(x,p ) = F^(m) = J«(x,z)[a 0 sin(/?z) - /?cos(/?z)]dz , (6a) 



, \ i/ ^ 2 °r / \ a 0 sin(Dz)- pcosipz) , x „ 

u(x,z) - F. \u) - — I u\x,p) 2 2 ~dp + S(x)e ** 

a o+P 



( 6 b) 



The subscript, s, stands for the smooth surface case and the last term in (6b), which 
represents a surface wave, decays with range and height and can be ignored for frequencies 
over 10 MHz according to Kuttler and Dockery [3], which will be the case for our work 
here. To better illustrate the fact this formula represents plane waves which are traveling 
towards and away from the boundary, it is rewritten in the mathematically equivalent form 
shown below: 



oo 1 °° 

u(x,p ) = F s (u) = ju(x,z)e ip: dz + — fu(x,z)e~ ipz dz } 

0 r s [p) 0 



(7a) 
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(7b) 



2 °° 

u(x,z) = F,(m) = - \u[x,p)\e~ ipz + r s (p)e ipz ]dp , 



where, r v (/?) 



P + fo o 
p-ia 0 



,p> 0, represents the plane wave reflection coefficient for smooth 



earth. From examination of (7b) it is apparent that the field is comprised of a wave traveling 
towards the boundary at z = 0, given by u(x,p)e~' p: , and one away from the boundary at z = 

0, given by T s (p)u(x,p)e ipz . It is evident that the reflected wave is equal to the incident 
wave multiplied by the reflection coefficient, which makes intuitive sense. 

Now that Fourier operators have been defined for the split-step algorithm in the case 
of a smooth ocean surface, we want to extend our derivation to include a boundary which 
consists of a rough ocean surface. Previous investigation of the problem of plane wave 
reflection from rough surfaces has resulted in the Kirchoff approximation and the concept 
of a rough surface reduction factor as given by Beckmann and Spizzichino [5]. To 
accommodate surface roughness for our problem we simply modify the smooth surface 
reflection coefficient, T s (p), by multiplying it with a rough surface reduction factor (defined 
below) to yield a rough surface reflection coefficient, T r (p). 

r r(p) = Po(A^)r i ( j p) . (8) 

The term, is the rough surface reduction factor and has been derived by Miller [2] 

as: 

Po{P’°h) = e ~ 2pV %( 2 P 2(y2 h) ■ ( 9 > 

The I 0 (-) term in (9) is the modified Bessel function of the first kind of order zero and a h 
represents the rms wave height deviation as determined by wind speed according to the 
formula: 
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c y h = 0.005 lfi 2 , 



( 10 ) 



where p = wind speed (m/s). A simpler form for the rough surface reduction factor has 
been recommended in CCIR report 1008-1 [6]: 




1 

V12Z-2 + V(32Z) 2 -7X + 9 ’ 



( 11 ) 



where X = 2p 2 a^ . 

We have now presented the necessary material to allow us to modify the Fourier 
operators such that the split-step algorithm can be used to predict propagation of vertically 
polarized waves over rough oceans surfaces. The field equation given for smooth surfaces at 
(7b) now becomes as follows for rough surfaces: 

1 °° 

u(x,z) = F;'[w(x, /?)] = — \u(x,p)[e- ipz + T r e ip: ] dp , (12) 



where, F~', represents the inverse Fourier operator for rough surfaces. The forward 
transform formula requires a bit more work to develop then simply substituting T r into (7a) 
for T s (as was done for the inverse transform) as this would not satisfy the consistency 

requirement f[f~‘(w)] = it . To derive the forward transform, both sides of (12) are 

alternately multiplied by e ±,pz , integrated with respect to z, and then appropriate linear 
combinations are taken. The result of this process yields the following for the forward 
Fourier transform: 

u(x,p)= Ju(x,z)e' pz dz + — 7 — r Ju(x,z)e~ ipz dz + i jk(p,q)u(x,q)dq , (13) 

o r r (p) o o 



where the kemal k[p,q ) in the last integral on the right hand side is: 



k(p,q) = 



r r (/?)-r r (g) 1 i-r r (p)r r (g) ' 

q-p q+p 



’27cr \p) * 



(14) 
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For p = q, a limit can be performed to yield: 





i -iM 



2 P 



1 

' 2*r ,{ P ) • 



( 15 ) 



where r/(p) represents the derivative of the rough surface reflection coefficient. In 
deriving (13), we have made use of the following identity from Papoulis [7]: 



je"** = n&(p) + - ! - , ( 16 ) 

o y 

where 5(-) is the delta function. 

Equation (13) may be rewritten in operator form by denoting the last integral on the 
right hand side as, K c u which yields: 



u(x,p) = (l-iK c ) 1 
= ¥ r [u(x,z)] 




( 17 ) 



where F r represents the forward Fourier transform for rough surfaces and / is the identity 
operator. For notational convenience we define the quantity P: 

P = {l-iK c ). ( 18 ) 
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III. SOLUTION PROCEDURE 



In the previous chapter we presented the case for using the parabolic equation and 
split-step algorithm to tackle the propagation problem of vertically polarized waves over 
rough ocean surfaces. We then defined a consistent forward and inverse Fourier transform 
pair to be used with the split-step algorithm. In this chapter we explain the solution 
procedure using this transform pair to determine field strengths at a given range. 

The solution procedure begins by recalling that we defined our transmitter source to 
be an omnidirectional point source at a height z = z v Numerically this initial field, «(x,z) at 
range x = 0, can be represented by a delta function. To determine the field at our first step 

downrange, (x = nAx) for« = 1, (17) is solved for u(x,p) with u(x,z) equal to the delta 
function. For this initial case, (17) becomes: 



u(x,p) = (l-iK c ) ‘ 



Jpz t 



7 ~*P*t 



r,M 



(19) 



Once the initial u(x,p) is known the u at the new range is determined by propagating in 

-ip 1 — 

free space. This is represented by the term e 2 *° in (2). Equation (12) can then be solved 
for the field u(x,z) at the first step downrange. A second phase correction is then applied to 
this field to compensate for variations in the index of refraction. This phase correction is 

seen as the term in (2). While this phase correction is suitable for most 

environmental scenarios, it is not considered adequate for rough surfaces when severe 
refractive index gradients are present. For this case we include the first and second 
derivatives of the refractive index and then modify (2) as shown below. The first and second 
derivatives were determined by doing a cubic spline interpolation of the given refractivity 
index data. The new formula that includes higher derivatives of the refractive index is: 
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ikgmAxi (m-l) 



n'Axf 



(m"Ax) 2 


0 Ax 




4 F"‘< 

r 


» ^ ■ 
k 

*9. 

1 


rs 




w 





( 20 ) 



where Ax = x - x 0 . Once a solution to (20) has been determined we can now take another 
step downrange, (x = nAx) for n = 2 , and repeat the process described above. This range 
stepping process continues until the field is determined at the desired distance downrange. 

For a given u{x,z) (17) is solved by numerically evaluating the quantity in the 
square brackets and then applying the inverse operator F [ on it. The numerical solution was 
obtained using MATLAB code to compute the Fourier transforms by means of an vV point 
complex FFT. In order to program the solution to (17) the continuous operator P first had to 
be put into a discrete form. This was done by using Simpson’s rule with weights, S n , to 
discretize the last integral in equation (13) to yield: 

a 

K c u(x,mAp) ~ A S„u(x, nAp)k(mAp, nAp) 

= Ku 



' N N' 

where K is of order [ — x — J and is the discrete version of the continuous K c . The upper 



N 

limit in the summation is — and not N because the integral is over the semi-infinite 



interval [0,°o) , whereas the complex FFTs assume limits of (- 00 , 00 ) . The lower limit is 1 
because u(x, 0) = 0 . The elements of K are: 

K mn = ApS n k{mAp, nAp ) , m, n = 1,2, . . . , ■ y , (22) 

As mentioned above, an vV-point FFT is used for the computation of the Fourier 
transform pair. Let us assume that the various quantities are band limited over -p mm < p < 
p max , and that the transform is evaluated ztp = 0, A p, 2 A p, . . ., (tV-1) A p, where: 
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( 23 ) 



4 P = 



2 % 

NAz 



Positive wavenumbers occur at p = while negative wavenumbers 

(N \ (N \ . N 

occur at + 1 — + 2 jAp,...\N -\)Ap . The value — Ap corresponds to both -p n 

and /W 



P 




To contain the computational domain vertically in the physical space, up to a 
maximum height z max , and to bandlimit the signal in the p-space, a Hanning window is 
employed in both the spatial and wavenumbers domains as follows: 
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h{n) = 



1 , 



sin 



Ann 

~JT'- 



A 3iV 

0<n< — 

3 N N 

——<n< — 
8 2 



Note that the Hanning window forces a gradual rolloff to zero of the data set as can be seen 
in Figure 4. 
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Figure 4. Hanning Window 
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IV. RESULTS 



In Chapters I - III we presented theoretical background and described the solution 
procedure for using the parabolic equation and Fourier split-step algorithm to solve the 
problem of propagation of vertically polarized waves over rough ocean surfaces. In this 
chapter we present numerical results for a variety of propagation scenarios. We will 
investigate the effects on propagation by varying parameters such as wind speed, step size, 
frequency, height of the upper boundary and refractivity index profile. Throughout this 
chapter we present plots showing Propagation Factor (PF), in dB, versus receiver height or 
distance downrange, where PF is defined as excess signal over free space. For a point 
source, the propagation factor can be computed from: 



where x is the downrange or horizontal distance from the transmitter, u is the field, and X 0 is 
the operating wavelength of the propagating signal. A positive (negative) value of 
propagation factor indicates a gain (loss) with respect to propagation in free space. 

A. CASE 1 - STANDARD ATMOSPHERE 

We begin our numerical results by running a test case against known data to validate 
our PE formulation. We look at propagation in standard atmosphere where the refractivity 

M = (m - 1) x 1 0 6 is given by: 



where z is the height in meters. The transmitter is at a height of z, = 30m, horizontal step 
size is Ax = 200m and the ground constants are given by e r = 80 and a = 4 Seimens per 
meter (S/m). Figure 5 shows the refractivity profile for the case of standard atmosphere. 
Figure 6 shows Propagation Factor versus Receiver Height at a range of 40 km with an 
upper boundary Z max = 512m, frequency of operation = 3 GHz and N= 512. This lobing 
pattern compares very favorably with that given by Kuttler and Dockery [3], Figures (7-10) 




(24) 



M{z) = (340 + 0.1 1 8z) , 



(25) 
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illustrate the effects of increasing wind speed and hence increasing surface roughness, for 
various heights of the upper boundary. Note that while the value of the upper boundary 
changes between Figures 7, 8 (150m) and Figures 9, 10 (300m), the vertical step size Az 
remains the same (0.146m). This is so because when the height of the upper boundary was 
doubled the number of points, N, in our FFT also doubled. For these figures, frequency = 
10 GHz, transmitter and receiver are both at a height of 25m, and Az = 0.146m. From 
Figures 7 and 8 one can see that as wind speed increases the excursions of the PF in the 
interference region is reduced and at longer ranges the PF curve no longer decays smoothly. 
The departure of propagation factor from monotonous decay for large ranges is due to 
numerical reflections from the upper boundary. Figures 9 and 10 have an upper boundary of 
300m. We can see by increasing Z max (and keeping the step size Az constant) the numerical 
reflections from the upper boundary do not affect the decay of the PF until further 
downrange. 



B. CASE 2 - TRI-LINEAR DUCT 



Our second case involves examining propagation where the index of refraction is 
characterized by a tri-linear duct as shown in Figure 11. Numerically, refractivity is given 
by: 



340 + 0.118z 0<z<135, 

M(z) = j 499.03 1.06z 135<z<150, 

322.33 + 0.1 18z z>150, 

where z is in meters. In this example the frequency of operation is 3 GHz, wind speed Ws = 
0, transmitter height = 30m, and Z max = 512m. We choose N- 512, resulting in a vertical 
step size Az = 2m, and choose a horizontal step size Ax = 200m. Figure 12 shows the 
Propagation Factor versus Receiver Height, at a range of 40 km. As for the case of a 
standard atmosphere, the propagation characteristics for the tri-linear duct are in agreement 
with Kuttler and Dockery [3]. 
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c. 



CASE 3 - EVAPORATION DUCT 



The next refractivity profile we investigate is that of an evaporation duct as depicted 
in Figure 13. The numerical values for the evaporation duct are shown in Table 1. 

Table 1. Refractivity Data For Evaporation Duct 



Height (m) 


M(z) 


0.000 


340.00 


0.135 


323.00 


0.223 


321.76 


0.368 


320.53 


0.607 


319.31 


1.000 


318.11 


1.649 


316.94 


2.718 


315.83 


4.482 


314.80 


7.389 


313.91 


12.182 


313.26 


20.000 


312.99 


33.115 


313.38 


54.598 


314.81 


90.017 


317.99 


148.413 


324.04 


165.000 


325.76 


300.000 


339.745 



Figures 14 and 15 illustrate that by wisely choosing the value of the upper 
boundary, accurate data can be obtained where it may not have been otherwise possible. 
Figure 14 shows propagation in an evaporation duct with Z max =75m, N=512, transmitter and 
receiver height = 25m, and ffequency=10 GHz. At the longer ranges the effects of the upper 
boundary and the evaporation duct prevent the PF from decaying in a stable manner. By 
increasing the upper boundary to 150m and N=1024 and keeping all others parameters the 
same. Figure 15 shows a stable PF can be achieved even at the longer ranges. 



D. EVAPORATION DUCT WITH ROUGH SEA SURFACE 

In this example we consider the same evaporation duct as above, but now we add 
the effects of a rough ocean surface. The wind speed considered is 10 m/s which results in a 
wave height of 0.51m from (11). Figure 16 shows both the propagation factors for a signal 
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in the evaporation duct with no wind Ws = 0 m/s and for Ws = 1 0 m/s. Other data are Z max = 
150m, frequency = 10 GHz and N = 1024. The effects of a rough ocean surface are more 
pronounced than the case of standard atmosphere. Specifically we see, (i) a reduction of the 
specular component as seen by the smaller excursions of the rough sea PF in the 
interference region, and (ii) increased losses for the rough sea PF over the smooth sea PF at 
greater distances. Figure 17 gives the results for PF vs. Receiver Height for the evaporation 
duct at a range of 100 km for wind speeds of 0 m/s, 10 m/s and 20 m/s. The upper 
boundary is Z max = 300m, frequency of operation =10 GHz, N=2048, and transmitter height 
z, = 25m. This figure illustrates the overall reduction in PF as surface roughness increases, 
however it also shows significant variation in the value of PF with respect to height due to 
the effect of the evaporation duct. 



E. 



CASE 4 - SURFACE DUCT 



The last example to be considered is the case of propagation in a surface duct. This 
refractivity is depicted in Figure 1 8 and its numerical form is given by : 

f350-0.335z 0 < z < 45.7 



M(z) = 



329.36 + 0.1 164z z>45.7 



Figure 19 shows the Propagation Factor vs. Range out to a range of 100 km for a signal 
traveling in the surface duct over smooth and rough ocean surfaces (Ws = 0 m/s and Ws = 
10 m/s). The frequency of operation is 10 GHz, transmitter and receiver height 
z, = 25m = z r , and an upper boundary Z max = 150m. As in the previous example, Figure 19 
illustrates that the rough ocean surface generally decreases the propagation factor as wind 
speed increases. The last example considered is Propagation Factor vs Receiver Height at a 
range of 100km for the surface duct with an upper boundary Z max = 300m. The results are 
shown in Figure 20 with all parameters remaining the same as the previous case except here 
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N = 2048. Once again surface roughness has reduced the propagation factor, and the effect 
of the surface duct has resulted in variations of the PF as receiver height is changed. 
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Figure 5. Refractivity for Standard Atmosphere. 
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PF(dB) 



Figure 6. PF vs. Receiver Height for standard atmosphere at range of 40 km. 

Transmitter height z, = 30m, wind speed Ws = 0, frequency = 3 GHz, N = 
512, vertical polarization and omnidirectional antenna. 
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Figure 7. Smooth and rough sea (Ws = 10 m/s) PF vs. Range for standard 

atmosphere out to a range of 100 km with upper boundary Z max — 150m. 
Transmitter height z t = 25m, receiver height z r = 25m, frequency = 10 GHz, 
N = 1024, vertical polarization and omnidirectional antenna. 
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Figure 8. Smooth and rough sea (Ws = 20 m/s) PF vs. Range for standard 

atmosphere out to a range of 100 km with upper boundary Z max = 150m. 
Transmitter height z, = 25m, receiver height z r = 25m, frequency = 10 GHz, 
N = 1024, vertical polarization and omnidirectional antenna. 
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Figure 9. Smooth and rough sea (Ws = 10 m/s) PF vs. Range for standard 

atmosphere out to a range of 100 km with upper boundary Z raax = 300m. 
Transmitter height z t = 25m, receiver height z r = 25m, frequency = 10 GHz, 
N = 2048, vertical polarization and omnidirectional antenna. 
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Figure 10. Smooth and rough sea (Ws = 20 m/s) PF vs. Range for standard 

atmosphere out to a range of 100 km with upper boundary Z max = 300m. 
Transmitter height z, = 25m, receiver height z r = 25m, frequency = 10 GHz, 
N = 2048, vertical polarization and omnidirectional antenna. 
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Figure 11. Refractivity for Tri-Linear Duct 



28 




PF(dB) 



Figure 12. PF vs. Receiver Height for tri-linear duct at range of 40 km. Transmitter 
height z t = 30m, wind speed Ws = 0, frequency = 3 GHz, N = 512, vertical 
polarization and omnidirectional antenna. 
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Figure 13. Refractivity for Evaporation Duct. 
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Range (km) 

Figure 14. PF vs. Range for evaporation duct to a range of 100 km with upper 

boundary Z max = 75m. Transmitter height z t = 25m, receiver height z r = 
25m, frequency = 10 GHz, wind speed Ws = 0, N = 512, vertical 
polarization and omnidirectional antenna. 
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Figure 15. PF vs. Range for evaporation duct to a range of 100 km with upper 

boundary Z max = 150m. Transmitter height z, = 25m, receiver height z r = 
25m, frequency = 10 GHz, wind speed Ws = 0, N = 1024, vertical 
polarization and omnidirectional antenna. 
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N = 1024, vertical polarization and omnidirectional antenna. 
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Figure 17. PF vs. Range for evaporation duct to a range of 100 km with wind speed = 
0 m/s, 10 m/s, 20 m/s. Upper boundary Z mltx = 300m, transmitter height z t = 
25m, receiver height z r = 25m, frequency = 10 GHz, N = 2048, vertical 
polarization and omnidirectional antenna. 



34 




Figure 18. Refractivity for Surface Duct. 
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Figure 19. Smooth and rough sea (WS = 10 m/s) PF vs. Range for surface duct out to a 
range of 100 km with upper boundary Z miu = 150m. Transmitter height z, = 
25m, receiver height z r = 25m, frequency = 10 GHz, wind speed Ws = 0 
m/s, N = 1024, vertical polarization and omnidirectional antenna. 



36 
















I 






PF (dB) 



Figure 20. PF vs. Receiver Height for evaporation duct at a range of 100 km with 

wind speeds Ws = 0 m/s, 10 m/s and 20 m/s. Upper boundary Z max = 300m. 
Transmitter height z t = 25m, receiver height z r = 25m, frequency =10 GHz, 
N = 2048, vertical polarization and omnidirectional antenna. 
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V. CONCLUSIONS 



In this thesis, a method of predicting radiowave propagation of vertically polarized 
waves over rough ocean surfaces was implemented and tested. An efficient solution was 
made possible by only considering one way propagation, thus allowing us to use the split- 
step parabolic equation method to tackle our problem instead of the much more 
computationally intensive Helmholtz equation. Our method incorporates the contributions 
from direct and reflected waves at any distance downrange from the transmitter, but ignores 
the effects of backscattering. This allowed us to develop a range stepping technique to 
determine the strength of the propagating waves. The fact that the split-step technique 
involves Fourier kernels (plane waves) made it possible to incorporate the effects of sea 
surface roughness directly into the spectral domain through the use of a rough surface 
reduction factor. 

The main purpose of this thesis was to modify the equations developed by 
Janaswamy [1] so the propagation of vertically polarized waves could also be solved by the 
split-step PE method. Numerical results for the case of vertical polarization show great 
similarity with those described in that report for the various refractivity profiles examined. 
The effects of the upper boundary layer are a significant factor in the accuracy of our model 
and must be taken into consideration when interpreting data. The user must take into 
consideration the value of the upper boundary, surface roughness and vertical step size as 
determined by the number of points, N, used in the FFT, to ensure accurate results are 
obtained with this model. 

The model presented in this thesis is applicable to propagation predictions for areas 
such as communications and radar. It is a useful tool for those designing, developing, and 
deploying new systems and for those that may be analyzing systems already in use. 
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